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1 Abstract 

In the flutter or forced response analysis of a 
turbomachine blade row, the blade row in question is 
commonly treated as if it is isolated from the neighboring 
blade rows. Disturbances created by vibrating blades are 
then free to propagate away from this blade row without 
being disturbed. In reality, neighboring blade rows will 
reflect some portion of this wave energy back toward the 
vibrating blades, causing additional unsteady forces on 
them. It is of fundamental importance to determine whether 
or not these reflected waves can have a significant effect on 
the aeroelastic stability or forced response of a blade row. 
Therefore, a procedure to calculate intra-blade-row unsteady 
aerodynamic interactions has been developed which relies 
upon results available from isolated blade row unsteady 
aerodynamic analyses. In addition, an unsteady 
aerodynamic influence coefficient technique is used to 
obtain a model for the vibratory response in which the 
neighboring blade rows are also flexible. The flutter 
analysis shows that interaction effects can be destabilizing, 
and the forced response analysis shows that interaction 
effects can result in a significant increase in the resonant 
response of a blade row, __ 

2 Nomenclature 

a speed of sound 

A nondimensional acoustic wave coefficient 

c chord 

G nondimensional gust coefficient 

i 

Im imaginary part 

M relative Mach number 

R nondimensional reflection coefficient 

Re real part 

s blade tangential spacing 

T nondimensional transmission coefficient 
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U rotor wheel speed 

V fluid velocity 

a angle of blade oscillation 

P inteiblade phase angle 

5 axial gap between neighboring blade rows 

r| tangential coordinate 

4 axial coordinate 

p fluid density 

© angular frequency 

© reduced frequency, coc/K 

3 Introduction 

In the flutter analysis of a turbomachine blade row, the 
blade row is commonly assumed to be isolated - 
disturbances created by the vibrating blades are free to 
propagate away from the blade row without being disturbed. 
Thus any reflections of these outgoing waves by other 
structural members or nonuniformities in the mean flow 
field are neglected. Although the forced response problem 
is by definition concerned with blade row interaction (with 
response to inlet distortion an important exception), forced 
response analyses also typically neglect any reflections of 
outgoing waves. However, in an engine environment, 
structural elements such as neighboring blade rows or struts 
and nonuniformities in the mean flow field will generally 
reflect some of this wave energy back toward the vibrating 
blades, causing additional unsteady forces on them. 
Whether or not these reflected waves can significantly affect 
the aeroelastic stability or forced response of a blade row is a 
question of fundamental importance. 

Several investigations have focused on the unsteady 
aerodynamic interaction between two rigid blade rows. Kaji 
and Okazaki (1970) investigated the unsteady aerodynamic 
interaction of two blade rows for the purpose of predicting 
rotor-stator interaction noise. They obtained a simultaneous 
solution to the unsteady lift distributions on both of the 
blade rows. Hanson (1992) modified Smith’s (1972) 
isolated blade row unsteady aerodynamic analysis to predict 
rotor-stator interaction noise, essentially extending Kaji and 
Okazaki's work to include effects of frequency scattering 
and mean flow turning by the blade rows. For a 
counter-rotating propfan in incompressible flow, Chen and 
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Williams (1991) used a panel method to determine the 
unsteady loads on rigid propeller blades. From the 
point-of-view of the present investigation, all of these 
investigations are somewhat limited because they did not 
delve into the aeroelastic problem and they were limited to 
two blade rows. 

One investigation that did consider the aeroelastic 
effects of interactions of a vibrating blade row and an 
adjacent structure was that of Williams, Cho and Dalton 
(1991). A three-dimensional linearized compressible panel 
method was used to calculate the unsteady aerodynamics of 
a ducted fan. The aeroelastic analysis allowed flexibility of 
both the fan and the duct. The duct was found to have a 
destabilizing effect on the fan. 

. s 

A number of time-accurate solutions to the Euler and 
Navier-Stokes equations for rotor-stator interaction have 
been presented, e.g., Rai (1987), Jorgensen and Chima 
(1989), Giles (1990). Due to the large computing time 
required to obtain solutions for just two blade rows, 
time-marching codes such as these are not practical for 
routine aeroelastic analysis purposes. 

The purpose of this paper is to investigate the effects of 
blade row interaction on flutter and forced response. A 
model for interaction between any number of blade rows is 
proposed which relies upon unsteady aerodynamic 
coefficients obtained from isolated blade row unsteady 
aerodynamic analyses. This allows the use of 
computationally-efficient linearized analyses, which are 
currently favored for turbomachinery aeroelastic predictions, 
for the analysis of multiple blade row unsteady aerodynamic 
interactions. 

The linearized analyses are generally capable of 
providing solutions for the unsteady blade forces and the 
outgoing waves generated by two classes of disturbances: 
(1) blade motion and (2) incoming waves such as acoustic 
and vorticity waves. Thus the unsteady aerodynamic 
coefficients which quantify the wave reflection and 
transmission characteristics of the (rigid) blade rows, as 
required to aerodynamically couple the blade rows, are 
provided by the linearized unsteady aerodynamic analyses. 
In addition, the response coefficients due to blade motion 
allow the calculation of influence coefficients which express 
the unsteady forces on one blade row due to the motion of 
another. By combining the various influence coefficients 
with equations of motion for each of the blade rows, a 
system of equations may be obtained which models what 
will be referred to as the 'dynamic coupling' of flexible blade 
rows. 

In this paper, a model for unsteady aerodynamic 
interaction between any number of blade rows is developed. 
The coupled equations of motion for the flexible blade rows, 


which rely upon unsteady aerodynamic influence 
coefficients for the intra-blade row coupling, are derived. A 
system of linear equations is found for the unsteady 
aerodynamic influence coefficients. Whitehead's subsonic 
flat plate cascade analysis, LINSUB (Whitehead, 1987), is 
used to provide the unsteady aerodynamic input for the 
influence coefficient solution although this model is suitable 
for use with the predictions of more sophisticated isolated 
blade row analyses which account for nonuniform mean 
flow fields (e.g., Verdon (1990), Hall and Clark (1991), 
Fang and Atassi (1993)). Before applying this technique to 
aeroelastic problems, the predictions of the unsteady 
aerodynamic model are compared with those of Kaji and 
Okazaki for outgoing acoustic wave amplitudes due to 
rotor-stator interaction. Then the model is used to 
investigate the effects of blade row interaction on the flutter 
and forced response of systems of three blade rows. 

4 Analysis 

4.1 Equation of motion 

A single degree-of-freedom model for torsional 
oscillation of an isolated blade row will be used to account 
for the effect of blade row interaction on the vibratory 
response of a row of identical blades. From Chiang (1988), 
the equation of motion for the small amplitude torsional 
oscillation of a two-dimensional section of a reference blade 
may be expressed as 

I^ + (\+2ig)Ia>la=M (1) 

a = a(0 is the angle of oscillation, I is the mass moment of 
inertia about the elastic axis, g is the structural damping 
coefficient for torsion and <o« is the undamped natural 
frequency for this mode of oscillation. M is the unsteady 
aerodynamic moment per unit span about the elastic axis. 

Simple harmonic motion at frequency © and a fixed 
interblade phase angle is assumed. Let the oscillation angle 
for a reference blade be_ a(r) = ae ml and the unsteady 
aerodynamic moment be M(t) = Me ,(S>t . Substituting these 
into Equation 1 and differentiating gives 

[-/© 2 + (1 + 2/g)/© 2 ]a =M (2) 

For the purpose of calculating M , the unsteady 
aerodynamic disturbances are assumed to be small relative 
to an inviscid, subsonic mean flow, resulting in a linear 
problem for the first order unsteady disturbance field. In 
this investigation, the unsteady disturbances are limited to 
acoustic and vorticity waves. Because the unsteady 
disturbance field is linear, the unsteady aerodynamic 
moment M may be determined by the superposition of 
motion-dependent and motion-independent parts. 
Typically, these moments are calculated under the 
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assumption that the blade row whose dynamic response is 
being analyzed is isolated from any other blade rows. 
Consider, for example, the response of a stator embedded in 
a multistage tuibomachine to the wakes of an upstream 
rotor. The unsteady aerodynamic coefficients are calculated 
based on a wake profile specified at the stator leading edge 
locus and the tangential velocity of the rotor. As a result of 
the unsteady interaction of the wakes with the stator, waves 
are created which propagate away from the stator: acoustic 
waves propagate both upstream and downstream, and 
vorticity waves convect downstream. If the stator is treated 
as an isolated blade row, as shown schematically in Figure 
1, these waves are free to propagate away from the stator 
undisturbed. In reality, these waves will impact neighboring 
blade rows, Figure 2. Even if the neighboring blade rows 
are perfectly rigid, they will emit additional acoustic and 
vorticity waves in response to the waves from the stator, and 
these waves, in turn, may further load the stator, which then 
emits more waves in response, and so on. Through these 
waves the blade rows are aerodynamically coupled. The 
question is, then, does this aerodynamic coupling between 
the blade rows have a significant effect on the aeroelastic 
response of the stator? 

One may cany this a step further and consider the 
vibratory response of the neighboring blade rows to the 
waves initiated by the stator/wake interaction. By allowing 
all of the blades to be flexible, they will respond by vibrating 
and emitting additional waves which will depend upon their 
vibrations. These motion-generated waves will further load 
the stator and result in additional waves ad infinitum. By 
allowing all of the blade rows to be flexible, the blade rows 
become dynamically coupled. Thus a second related 
question is, can the dynamic coupling of the blade rows 
significantly affect the vibratory response of the primary 
row? 

Since the unsteady aerodynamic problem has been 
assumed to be governed by a system of linear equations, the 
aerodynamic blade row interactions can be expressed as 
nondimensional aerodynamic influence coefficients. Let 
Cjk be the nondimensional moment coefficient for the y-th 
blade row due to oscillation of the £-th blade row which 
includes effects due to wave interaction (i.e., reflection and 
transmission of acoustic waves, shed vorticity) with the 
other non-vibrating blade rows. Thus C* is the moment 
coefficient on the y-th blade row due to its own oscillations. 
These coefficients can be used to dynamically couple any 
number of blade rows. For a system of n blade rows 
restricted to torsional oscillation, Equation 2 will be applied 
to each blade row. The unsteady aerodynamic influence 
coefficients are used to express the moment in Equation 2 as 
the linear combination of the influences of all of the blade 
rows. Assuming that the forced response of the y-th blade 


row is being investigated, the unsteady aerodynamic 
moment on a reference blade of that row is 

Mj - pjV]c](Cj\CL] + CjiCLi + ... + Cjj a J + + C>,a fl + 

[CGlfi + [CG] j7 % + ... + [CG]4 + ... +[CG]^). (3) 

Pj is the gas density, Vj is the relative inlet velocity and c ; 
is the chord for blade row y. a* is the complex-valued angle 
of oscillation of the reference blade of row k. [CG] ;jt is the 
nondimensional moment coefficient for row j due to a gust 
acting on row k with upwash w k - typically [CG] ;ir is used to 
model the effects of viscous wakes. Substituting Equation 3 
into Equation 2 then rearranging into a nondimensional 
form results in 


Cyiai +C,20t2 + ... 
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TT 2 2 f 


... +C/n-l)(X»-l +C,na„ = -Z 7^-[CG] * (4) 

b= 1 V k 


- CO/C/ / h 

where 5, = is the reduced frequency, r, = j is 

J Vj y V 

the nondimensional radius of gyration of the blade section 

and j ij = — — — - is the mass ratio of that section. Letting 
7tp(c ; /2) 

where y j = Equation 4 becomes 

C/i<Xi + ... +(Sy H- Cj r )Oy+ ... +Cj n CL n = — L “’[CG] ; * (5) 


The only relevant values of [CGJ^are those with the 
frequency <ay . In many cases, the gusts acting at any row 
k*j will not result in excitation of row j at ©y, hence 
[CG] Jk = 0 for k * j and Equation 5 becomes 

Cyiai + ... + (Sy + Cy/)ay + ... + Cy n a n = — jr[CG]^ . (6a) 

In this case, the equation of motion for a blade row k *y is 
Cjti ai + ... h- (S jt Cjfc*)a.jt 4- ... 4- = — j^[ CG] kj . (6b) 

Using Eq. 6a for blade row j and Eq. 6b for the 
remaining n-1 blade rows, a linear system of equations is 
obtained which may be solved for the dynamics of these 
blade rows due to a gust acting on row y. This system may 
be expressed as 

Ax=B (7) 

where 
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x = [ai,a 2 . 

.. 

.,a^i,a n ] T 

(9) 

B = -[ 

g[CG]„, 

Wj 

V} 

[CG] 2y >.. 

.^[CG] Jp .... 



?[c<Viv^ CG W T < lo :> 

Specifying B=0 in Equation 7 results in a homogeneous 
system appropriate for determination of the system stability: 

Ax=0 (11) 

A meaningful solution exists only if the determinant of the 
coefficient matrix A is equal to zero. 


flutter. Im(/V JyJ) is the ratio of the blade row j oscillation 
frequency to its natural frequency. 

4.2 Coupled blade row unsteady aerodynamics 

Expressions for the unsteady aerodynamic influence 
coefficients will now be developed. The mean flow field is 
assumed to be subsonic, and only those acoustic waves 
which propagate unattenuated will be considered. Blade 
row geometry definitions and flow field definitions are 
shown in Figure 3. 


Assume that blade row j is oscillating harmonically 
with unit amplitude at frequency ©, and interblade phase 
angle )3y. Due to these oscillations, acoustic waves 
propagate both upstream and downstream, and vorticity 
waves are convected downstream. As shown in Appendix 
A, acoustic waves will propagate unattenuated through a 
uniform flow field of Mach number M<\ only if the 
tangential wavenumber m (i.e., the r\ component of the 
wavenumber) satisfies 


m 


< © 


> a 

1 -M 2 




J 


(14) 


For any blade row k, the structural coefficients 
r*, and gjt are assumed to be known. To investigate the 
stability of a particular row, row j 9 © j is chosen and this, 
along with the geometry and mean flow field, sets the values 
of 5* for the other blade rows (see Appendix B). The 
frequency ratios y* are unknown, but they may be expressed 
as a functions ofy,: 

7* 2 2 2 



( 12 ) 


and are the § and r\ components of the Mach 
number, © is the frequency and a is the speed of sound. If 
this condition is met, two propagating acoustic waves will 
be produced with tangential wavenumber m , one 
propagating upstream and the other downstream. The (-) 
sign in Equation 14 corresponds to the upstream-running 
wave and + corresponds to the downstream-running wave; 
the < symbol corresponds to the + and the > corresponds to 
the (-). These waves are said to be superresonant or cut-on. 
Otherwise, the waves decay exponentially with axial 
distance and are said to be subresonant or cut-off. To be 
consistent with the interblade phase angle (3 j9 m must satisfy 

mSj - p; + 2tct, r = 0, ±1, ±2,... (16) 


co u ic and ©„, are assumed to be known, and the ratios ©,/©* 
will also be known due to their dependence on known 
parameters. Thus the determinant of the coefficient matrix 
yields a polynomial of degree n in y ) which may be solved 

for Y;. 

The system stability is determined by the values of y 
Since 


and aj{/)=a Rc(i/Jyjy>0 indicates growing 
oscillations for blade row j and blade rows k (through the 
dependence of ©jt on ©, ) and therefore system instability or 


where 5) is the spacing of the row j blades. As shown in the 
Appendix, scattering of these waves at the other blade rows 
results in waves returning to row j which are generally not 
at the frequency of interest. Although these waves might 
eventually be scattered back to the frequency of interest, the 
scattered waves will be neglected. 

Let the waves generated by the motion of row j be 
denoted Gmj for the gust or vorticity wave, AJjj for the 
upstream-traveling acoustic wave and A^j for the 
downstream-traveling acoustic wave. Superscripts of (-) and 
(+) are used to indicate the directions of acoustic waves, (-) 
being a wave going toward -oo (upstream) and (+) being 
toward -too (downstream). To avoid unnecessary 
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complications in the present explanation, acoustic waves are 
assumed to propagate unattenuated in both the upstream and 
downstream directions for only one tangential wavenumber 
consistent with the interblade phase angle. Assuming that a 
dynamic steady state is reached for all the resulting waves, 
each blade row k will have the net effects of three different 
incoming waves acting upon it, a gust, G/,*, an 
upstream-running acoustic wave, and a 

downstream-running acoustic wave, A . The outgoing 

waves at row k are linearly related to the incoming waves by 
the various reflection and transmission coefficients of the 
blade row. For example, the outgoing upstream-running 
acoustic wave will be 

A o.k +A^ Jt k +G Jt kA G Jt (17) 

where, for row k, T k is the transmission coefficient for 
upstream-traveling acoustic waves, R k is the reflection 
coefficient for downstream-traveling acoustic waves and 
A Gk is the upstream-running acoustic wave due to an 
incoming gust Similarly, the other two outgoing waves are 
a linear combination of the incoming waves: 

A O t k = A I t k^k Tjc +Gi f icAQ k (18) 

and 

G 0 ,k -A~ } k G~ A k +-A£ k G A k +G/,*Gg,*. (19) 

The new coefficients G~ Ak , G\ k and Ga.^are the resulting 
gusts or vorticity waves due to an upstream-running acoustic 
wave, a downstream-running acoustic wave and a gust 
respectively, incident on row k. 

By defining the outgoing and incoming waves at 
common interfaces between the blade rows, the outgoing 
waves at row k are related to the incoming waves at adjacent 
blade rows by 

A 0,k ~ A J t k-i 

J+ — A + 

A 0,k Jt+1 

Go,k = Gi k + 1 

so that Equations 17 through 19 become 

A J,k-\ ^ A J P k^k +A J,k^k + Gj'j c A C j c (20) 

A J, Jt+1 -^I.lftk + A ],k Tk +G] t l c A G j c (21) 

.. . Gj t jt+i =Aj tk G Ak +AJ k G A k + G] tk G G ,.k (22) 

For the blade rows neighboring the oscillating row, the 
incoming waves may be split into known parts due to the 
oscillation and unknown parts, indicated by hats, due- to the 
ensuing interactions with the other blade rows: 


A 

Alj - 1 = ^/.j-l + ^Mj 
^ Jj+l = -d/j+i 

A 

Gi t j + 1 = Gi y p. i + Gmj 

Then, for row j and its immediate neighbors, the system of 
equations is, with all of the unknowns collected on the right 
hand side and the known terms due to the motion of row j 
on the left hand side, 

a mj = +A i.jTj +A]jRj + GijA G j 

0 = -A*JJ_ X +A]j_ 2 Rj-2 +>*/,^ 2^2 

0 = -Gl t j - 1 +Ajj_ 2 G A j. 2 2^A,J-2 + G/,/~2Gg,J-2 

i -Gm,jAq^ } =-Alj+A] J+] T~+ ] + 

\Rj+\ + GlJ+\A QJ + 1 

=-Al J +AJj_ l R~_ ] +A;j_ l T*L l +Gj t j-\A G ^ (23) 

A 

-d MJ G A 'j _ j = -Gj t j +Aj j_ ] G a +A 1 ^G A j_< [ + 

Gi t j-\Gc,j- i 

0 = 2^>f2 + A lj+ 2^1+2 + Gl t j+2 A C j+2 

A Mj = ~ A lj+\ + A i,jRj + Gi'jA G j 

Ghj = -Gi t j+\ +A]jG A j +Aj tJ G A j + GjjGgj 

Assuming that three blade rows will be analyzed, the 
present system of 9 equations has 15 unknowns because the 
incoming waves to rows (/- 2) and (/+ 2) are required. By 
assuming that no waves internal to the system are reflected 
by these two rows, = G Aj _ 2 =Rj+ 2 =a g,j + 2 = °> and 
that there are no waves incoming from sources external to 
the system, Ajj_ 2 = Gj.j- 2 = A ij +2 = °> the number of 
unknowns is reduced to 9. Thus the system of equations 
may be solved for the incoming waves to each of the three 
blade rows. Using coefficients provided by the isolated 
blade row analysis which relate the moment on a blade to 
the incoming waves, the unsteady moment on each blade 
row k due to oscillation of row j - the influence coefficient 
Ckj - may then be determined. Similarly, by using the 
waves generated by a unit amplitude gust acting on row j as 
the known terms in Equations 23, the gust coefficients 
[CG]* ; may also be determined. 

In this investigation, Whitehead's (1987) linearized 
subsonic flat plate cascade code, LINSUB, which is based on 
Smith (1972), will be used to supply the required 
aerodynamic quantities. It is based on assumptions of 
inviscid, isentropic, subsonic flow through an isolated, 


6 


5.2 Flutter 


infinite cascade of equally-spaced flat plate airfoils. The 
airfoils are set at zero mean incidence so that the mean flow 
is uniform, and the unsteady disturbances created by the 
airfoil oscillations or incoming gusts are assumed to create 
small unsteady disturbances to the mean flow, resulting in a 
linear problem for the unsteady flow. LINSUB calculates 
moment coefficients due to oscillation of the cascade or due 
to incoming acoustic or vorticity gusts, and outgoing 
acoustic and vorticity waves due to oscillation of the cascade 
or due to incoming waves. Thus all the coefficients required 
by Equation 23 can be obtained from LINSUB. 

5 Results 

This model for unsteady aerodynamic blade row 
interaction will be applied to both the flutter and forced 
response of the middle blade row in a system of three blade 
rows. Before investigating these aeroelastic problems, the 
predictions of the unsteady aerodynamic model will be 
compared with the predictions of Kaji and Okazaki (1970) 
for the outgoing acoustic wave amplitudes due to 
rotor-stator interaction. 

5.1 Acoustic wave generation by rotor-stator 
interaction 

Kaji and Okazaki investigated the unsteady 
aerodynamic interaction of two rigid blade rows in motion 
relative to one another for the purpose of predicting 
rotor-stator interaction noise. The mean flow was a 
uniform, compressible stream about two cascades of flat 
plates aligned with the flow. They obtained a simultaneous 
solution to the unsteady lift distributions on both of the 
blade rows and determined the resulting outgoing acoustic 
waves. For excitation of the downstream blade row by 
wakes from the upstream row, the wake model of 
Silverstein, Katzoff and Bullivant (1939) was used to specify 
the wake velocity distribution at the leading edge of the 
downstrea m b lade row. 

The stagger angles of the blade rows are 
yj = — y 2 =30°, giving identical relative Mach numbers 
for the two rows. The solidity of each row is 1 and the 
reduced frequency is 7C. Comparison between the 
predictions of the present method and those of Kaji and 
Okazaki is shown in Figure 4. The amplitudes of the first 
harmonic upstream and downstream acoustic waves are 
shown as a function of the relative Mach number into the 
blade rows and the axial gap between the blade rows as a 
fraction of chord. There is generally good agreement 
between the predictions for both of the axial gaps. 


To investigate the effects of unsteady aerodynamic 
blade row interaction on flutter, two configurations of three 
blade rows were analyzed. The first configuration consists 
of a staggered rotor placed between unstaggered stators 
while the second has an unstaggered stator between 
staggered rotors. To help isolate various effects, the stability 
of the middle blade row will be determined for three 
different cases: (1) the middle row is treated as if it is an 
isolated blade row, so there is no unsteady aerodynamic 
coupling between the blade rows (the NC' solution); (2) 
unsteady aerodynamic coupling occurs between the middle 
row and its neighbors, but the neighboring blade rows are 
rigid (the 'AC solution); (3) unsteady aerodynamic 
coupling occurs and the secondary blade rows are now 
flexible, hence the blade rows are dynamically coupled (the 
VC solution). Variation of the blade row spacing as a 
technique for enhancing blade row stability will be 
investigated. 

System parameters for configuration I are shown in 
Table 1. The stators are unstaggered and the axial Mach 
number relative to the stators is 0.4. Because the relative 
flow is required to be aligned with the rotor and the rotor is 
staggered 60°, the rotor relative Mach number is 0.8. Rotor 
velocity triangles were then used to find the rotor blade 
Mach number, £//<z=- 0.693. For. the rotor geometry, mean 
flow conditions and reduced frequency 5 2 = 1 , the acoustic 
resonances occur at = -20.3° and P 2 - 146,4° . The radii 
of gyration and mass ratios are similar to those used by 
Chiang (1988), and are believed to be reasonable values for 
titanium compressor blades. For the initial calculations, the 
leading edge loci of the blade rows are at axial locations 
-1.0, 0.1 and 1.5 chord, so there is a 0.1 chord axial gap 
between the rotor and the upstream stator and a 0.9 chord 
axial gap between the rotor and the downstream stator. 
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Table 1 Blade row configuration I 



Stator 1 

Rotor 

Stator 2 

Blade row index, h ill 

1 

2 

3 

Chord,*?"'" 1 

1 

1 

1 

Solidify, M 

1.5 

1.4 

1.5 

Stagger angle (deg.)1|||' 

0 

60 

0 

Relative Mach no., M 

0.4 

0.8 

0.4 

Blade Mach number, V/a 

0 

-0.693 

0 

Elastic axis / chord "Jiff- 

0.4 

0.4 

0.4 

Radius of gyration 

0.4 

0.4 

0.4 

Mass ratio 

200 

200 

200 

Natural freq., ©»/©» Rowe 

1 

1 

1 

Damping coefficient 

0.5% 

0.5% 

0.5% 


To obtain an equation of motion for the isolated rotor, 
the oscillations of the blade rows are decoupled by 
modifying the coefficient matrix A* In Equation 8, letting 
/= 2, n = 3 and C,* = 0 for all k * j then substituting this A 
matrix into Equation 1 1 results in 

(S 2 + Cz2)<X2 = 0 (24) 

for row 2. Thus the oscillations of the blade rows are 
decoupled, but Equation 24 still retains the effects of 
unsteady aerodynamic interaction between the blade rows 
through the coefficient C 22 • Replacing C 22 by the isolated 
blade row motion-induced moment, C 2 , results in an 


equation for the isolated rotor: 


(S 2 + C 2 )a 2 = 0. 

(25) 

y 2 is determined by solution of 


s 2 + c 2 = 0, 

(26) 

then the oscillation frequency 

is found from 


7g>2/©«2 = H JW . Solutions obtained for the isolated blade 
row response (NC solution) for 15° increments in (3 2 are 
shown in Figure 5; the corresponding value of (3 2 in degrees 
is given adjacent to each point. The rotor is stable for all p 2 
(Re(/a) 2 /© tt2 ) <0). Subresonant values of the interblade 
phase angle (p 2 <-20.3° and p 2 > 146.4°) are the most 
stable and form a relatively tight group of points near 
Re(/co 2 /co ui) = “0. 12 . Except for p 2 = 1 50° , the 

superresonant points form another group which is much 
closer to the neutral stability line. In particular, the points 
60° < p 2 £ 120° are relatively close to being unstable. 

When aerodynamic coupling is taken into account, but 
the neighboring blade rows are rigid (AC solution), 


Equation 24 shows that the rotor frequency-of-oscillation 
must satisfy 

S 2 + C 22 = 0. (27) 

However, for subresonant values of p 2 , all of the acoustic 
waves decay with axial distance and are therefore neglected 
in the analysis. Thus in the subresonant region the solutions 
to Equation 27 will not differ from those to Equation 26. 
The solutions for the superresonant interblade phase angles 
have changed, Figure 6, and some of the roots indicate 
flutter - the points 75° ^ p 2 < 135° have all shifted right far 
enough that Re(/© 2 /©« 2 ) > 0. 

Solutions to Equation 11 without any simplification 
(aerodynamic coupling with flexible neighboring blade 
rows, the DC solution) are shown in Figure 7. About 5 
minutes of computing time on a 80486 personal computer 
was required to obtain these solutions. Now that all of the 
blade rows are allowed to oscillate, there are three roots at 
each superresonant intefblade phase angle. Where two roots 
are difficult to distinguish from one another, *(2)' is placed 
next to the value of p 2 . Two distinct groups of roots are 
shown in the figure - one group, with Im(/© 2 /© M2 )a 1, 
stretches horizontally across the plot while the second group 
stretches vertically just to the stable side of the neutral 
stability line. The first group is replotted in Figure 8 along 
with the AC solution first shown in Figure 6. These two 
solutions are effectively identical, implying that the 
dynamics of the stator blades have a negligible effect on the 
rotor stability. Upon inspection of the coefficient matrix A 
for this configuration, it was found that the off-diagonal 
terms were small relative to the diagonal terms, thus the 
oscillations of the blade rows were effectively decoupled due 
to weak unsteady aerodynamic interaction. Due to the 
decoupling, the remaining solutions - the vertical group of 
roots in Figure 7 - appear to be AC solutions for the stability 
of the two stators. But the solutions for the stators are 
nearly identical despite their differing locations in the 
configuration, which indicates that unsteady aerodynamic 
coupling effects on them are also negligible. The 
frequencies of these solutions vary with interblade phase 
angle because they are presented in terms of the rotor 
frequency - the stators are oscillating near their natural 
frequencies, but due to the motion of the rotor, the frequency 
in the rotor frame varies with the interblade phase angle or 
tangential wavenumber (Appendix B). 

The rotor stability will vary with spacing between it and 
the upstream stator, 5. Sample results for the present 
configuration are shown in Figure 9; the positions of the 
stators are fixed while the rotor position between them is 
varied between 8/c=0 (rotor leading edge coincident with 
upstream stator trailing edge) and 5/c-l (rotor trailing edge 
coincident with downstream stator leading edge). Relative 
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to the isolated rotor solution, the unsteady aerodynamic 
coupling effects are generally destabilizing, but for 
5/c>0.68, a slight stabilizing effect occurs. The rotor is 
unstable for 5/c<0.31. These results suggest that blade row 
stability may be enhanced by variation of the blade row 
spacing. 

The stability of a second configuration was also 
analyzed; configuration n. Table 2, consists of an 
unstaggered stator between two identical rotor blade rows. 
The blade row spacing was chosen so that a 1/2 chord axial 
gap separates the stator from each of the neighboring rotors. 
For the stator geometry, mean flow conditions and reduced 
frequency ©2 = 1, the resonant interblade phase angles are 
±19.2°, so the superresonant region is relatively small. 
Figure 10 shows the NC solution for the stator along with 
the AC solution. The stator is stable in both cases although 
aerodynamic interaction makes it less stable for some p: . 
However, if the rotors are allowed to vibrate, there are 
dramatic changes in the solution. From Figure 11, the DC 
solution has two unstable values of P2, 15° and 10°, and 
introduces a number of other roots which are quite different 
from those for the rigid rotors solutions. The dynamics of 
the neighboring blade rows, which were negligible for 
configuration I, play a crucial role in the stability of 
configuration II. 


Table 2 Blade row configuration IT 


1 : §i| | 

Rotor 1 

Stator 

Rotor 2 

Bladerowindex,£ 

1 

2 

3 

Chord, ’ v • 

1 

1 

1 

Solidity, c/s 

1.5 

1.3 

1.5 

Stagger angle (deg.) 

60 

0 

60 

felatiye ]^ch h6. r A/ 

0.8 

0.4 

0.8 

Blade Mach no., U/a ,, 

-0.693 

0 

-0.693 

Elastic axis / chord 

0.5 

0.5 

0.5 

Radius of gyration 

0.4 

0.4 

0.4 

Mass ratio 

200 

200 

200 

Natural freq., ©u/©vsutor 

1 

1 

1 

Damping coefficient 

0.5% 

0.5% 

0.5% 


5.3 Forced response 

To investigate the effects of unsteady aerodynamic 
blade row interaction on forced response, calculations were 
made for the first harmonic response of Jhe stator in 
configuration II due to the wakes of the upstream blade row. 
The first harmonic of the wake is assumed to have unity 


amplitude. The stator reduced frequency is ©2=16.32. In 
contrast to the stability solution procedure, where ratio of 
the oscillation frequency to the undamped natural frequency 
was the solution, the oscillation frequency is now fixed and 
the undamped natural frequency will be varied in order to 
see how the response varies in the vicinity of resonance. 

The forced response of the configuration II stator is 
shown in Figure 12 as a function of ©2/©«2> the ratio of the 
stator excitation frequency to the stator undamped natural 
frequency. The response parameter is the ratio of the 
oscillation amplitude with unsteady aerodynamic coupling 
effects to the oscillation amplitude for the isolated stator. At 
resonance, the AC response is nearly 45% greater than the 
NC response. Away from resonance, the AC response is 
still about 20% greater than the NC response. For this 
configuration, the DC response was essentially identical to 
the AC response - the oscillations of the rotors did not affect 
the stator. 

A brief parametric study of the effects of blade row 
spacingancTrotor solidity on the stator forced response was 
made. Configuration II was used for this study along with 
two configurations identical to it except for changes in the 
rotor solidities. The gap between the stator and the 
upstream rotor was varied between 0 and 1 chord; with a 
gap of one chord, the stator trailing edge was coincident 
with the leading edge of the downstream rotor. As shown in 
Figure 13, the response of the stator can vary greatly with 
these parameters. With the higher solidity rotors ( 
(c/s) 1 = (c/s ) 3 = 1.55), the stator response amplitude varies 
between exceeding the isolated stator response by nearly 
60% to being 15% less than the isolated stator response. In 
contrast, with the lower solidity rotors ((c/s) ]= (c/s ) 3 
=1.45), the stator response varies little with axial gap, 
staying at about 13% greater than the isolated stator 
response. The configuration II solution falls between these 
two in both maximum response and fluctuation of the 
response with blade row spacing. Although it might be 
tempting to draw some conclusions from the trends of these 
results, the very limited nature of this study, not to mention 
the complexity of cascade unsteady aerodynamics, prohibit 
this. 

6 Summary and Conclusion 

The effects of blade row unsteady aerodynamic 
interaction on flutter and forced response have been 
investigated. A model for interaction between any number 
of blade rows was developed which relies upon unsteady 
aerodynamic coefficients obtained from isolated blade row 
analyses. Using influence coefficients which express the 
unsteady forces on one blade row due to the motion of 
another, a structural model was obtained which accounts for 
coupling of the vibratory responses of multiple blade rows. 
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or dynamic coupling between the blade rows. A special case 
of this model is obtained when neighbors of the blade row 
under consideration are assumed to be rigid, but the 
unsteady aerodynamic coupling between these neighbors is 
retained. 

This analysis technique was applied to two model 
configurations, each consisting of three blade rows. 
Stability solutions were found for the middle blade row of 
each configuration. For the first configuration, 
aerodynamic coupling had a destabilizing influence relative 
to the isolated blade row solution - the isolated blade row 
solution indicated stability but the unsteady aerodynamic 
coupling solution indicated flutter. The additional effect of 
dynamic coupling was negligible relative to the 
aerodynamic coupling effect. In contrast, the dynamic 
coupling analysis indicated flutter for the second 
configuration while the unsteady aerodynamic coupling 
solution and the isolated blade row solution both predicted 
stability. It is concluded that both unsteady aerodynamic 
coupling and dynamic coupling between blade rows can 
have significant effects on blade row stability. 

The forced response analysis was applied to one 
baseline configuration. Significant increases in the response 
amplitude were found to occur due to unsteady aerodynamic 
coupling. A parametric study showed that blade row 
spacing and solidity can have significant effects on the 
response amplitude. 
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Appendix A Linearized theory of cascade 
wave propagation 


The unsteady aerodynamic blade row interaction 
analysis presented in this paper depends upon the ability to 
predict the characteristics of waves traveling between the 
blade rows, i.e., the amplitudes, phases and wavenumbers of 
these disturbances must be known. While the wave 
amplitudes must be determined by a numerical solution for 
the isolated blade row unsteady aerodynamics (see 
Whitehead (1987), Verdon (1990)), which is beyond the 
scope of this discussion, the wavenumbers and hence wave 
properties such as spatial decay rates may be determined 
with some analysis. 

The flow is assumed to be inviscid, isentropic, 
two-dimensional and subsonic. In addition, the mean flow 
is assumed to be uniform; although this assumption may 
appear questionable, it is consistent with current linearized 
cascade analyses which assume that regions of uniform flow 
exist at some distance upstream and downstream of the 
cascade. Unsteady disturbances are assumed to be small 
perturbations to the mean flow field, resulting in a sy stem of 
first order linear partial differential equations for the 
unknown perturbation quantities. When harmonic solutions 
which satisfy cascade periodicity are assumed, a system of 
linear equations is obtained, which may be solved for the 
disturbance wavenumbers. 

The mass conservation equation is, for the 
coordinate system of Figure 3, 

^ + p fg + ^l =0 . 

Dt dr\) 

The momentum equations are 

Du__ \$P 
Dt 9 81 


(A-l) 


(A-2) 


and 


Dv 

Dt 


\8p_ 

'pan' 


(A-3) 


are substituted into the governing conservation equations. 
Once the equations are expanded and products of 
perturbations are neglected, a system of linear partial 
differential equations results with the perturbation q uantit ies 
as the dependent variables. Conservation of mass becomes 


Dq£\ 

Dt 




The momentum equations become 


and 


with the operator 


DoU\ 

1 dp i 

Dt 

po as 

Dov\ _ 

1 dpi 

Dt 

Po aq 


(A-4) 


(A-5) 


(A-6) 


Do _ 8 . 8 ^ a 

Dt 8t as aq 


To close the system of equations, the speed of sound for 
isentropic flow of a calorically perfect gas, a, is introduced. 


2 _(&py _<*P_ 

a {dpj s dp 


(A-7) 


Expanding Equations A-4 through A-6 and using A-7 to 
replace derivatives of pi with derivatives of p i gives 


where 


a +A a +B ^=o 

dt 8% dr | 


q = (pi,»i,vi) 


(A-8) 


A = 


and 


B = 


u is the S component of velocity, v is the n component of 
velocity, p is the density and 

D_ S £ +U _L +V ± 

Dt dt 8$ 8r\' 

Using the subscript 0 to denote the constant mean values 
and 1 for the perturbations, 

p = po + pi, U = U 0 + U\, ... 


i/o poa 2 0 ^ 
1/po u o 0 

0 0 i/o J 


v 0 0 po<r 

0 vo 0 
1/po 0 vo 


The dependent variables are assumed to depend 
harmonically on the spatial position and the time. Thus the 
pressure perturbation is expressed as 

p x =p7e'(“ ,+ ^ + ' m i) 

where p] is the pressure disturbance amplitude, i = J-\ 
and / and m are the axial and tangential wave numbers, 
respectively. Analogous relations are used for u i and vi . 
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Substituting the perturbation expressions into Equation 
A-8, a system of linear algebraic equations is obtained: 


© + UqI + Vq/W 

a J /p 0 

2 > 
a'mpQ 

( — A 
P i 

Up 0 

© + «o/ + V(j7n 


«7 

mfpo 


ffl + UqI + Vq/W j 



For nonzero pressure disturbances, Equation A- 13 
implies <o + wo/+vo/n*0 which is true of the second 
solution. Equation A-12. From Equations A-17, mu\ = lv\. 
Substituting this into Equation A-14, it follows that the 
vorticity is zero; the second solution therefore corresponds 
to irrotational pressure perturbations. 


To obtain a non-trivial solution to this linear system, the 
determinant of the coefficient matrix in Equation A-9 is set 
equal to 0: 

(© + uol + vo/n)[(co + uqI + vo m ) 2 - a 2 (l 2 +m 2 )] = 0. (A-10) 
The two solutions to Equation A-10, 

cd + uoZ + vq/w = 0 (A-ll) 


and 

(© + wo/ + vo/n) 2 = a 2 (l 2 + m 2 ) (A-12) 

encompass different physical phenomena. 

From the momentum equations, relations between 
pressure and velocity fluctuations are obtained: 

p\ = ~y-(© + uol + vo/n)wT 

= -^-(<0 + UqI + V 0 m)v7. (A-13) 


Hence for the solution of Equation A-10 obtained specified 
by Equation A-ll, the pressure fluctuations are, from 
Equation A-13, zero, but the vorticity is 


u, gvi du\ 


(A- 14) 


The phase velocity of this disturbance is determined by 
tracking a constant disturbance phase. In this reference 
frame, the phase does not change with time: 

^(co/+^ + /nri) = 0 


or 

® + + = 0. (A- 15) 

dt dt 

Combining Equations A- 15 and A-ll to eliminate co results 
in 

l^j- = uol + VQtn. . . (A-16) 

dt dt 


Equation A-12 may be solved for the axial wave 
number: 


uo(® + vom)±a ]((£> +vom ) 2 -(a 2 -ul)m 2 
1= y — . (A-17) 


a 2 -u o 


The tangential wavenumber must satisfy cascade periodicity, 
and therefore be compatible with the interblade phase angle 
|3. Hence 


ms = P + 2nr, r= 0, ±1, ±2,... 


(A- 18) 


where s is the cascade spacing. Because m depends upon 
real quantities, it is also real. 


The nature of acoustic wave propagation depends upon 
/. From Equation A-17, if (a> + vom) 2 -(a 2 -ul)m 2 < 0, 
then / will be complex. Let l = l Re + il Im so that 

P\ = 

(A- 19) 

which indicates that the wave either grows or decays 
exponentially with axial distance. The amplifying wave 
solution is commonly disallowed on the basis that it is not 
physically acceptable (Fang and Atassi (1993), Verdon 
(1987)), leaving only the decaying wave solution. If 
(© + vo/w) 2 — (a 2 - u\)m 2 > 0, / is real and the acoustic 
waves will propagate unattenuated, one upstream (for the (-) 
sign in Equation A-17) and one downstream (for the + sign 
in Equation A-17) - see Whitehead (1987). At the acoustic 
resonance or cutoff condition, (co + vo/n) 2 =(,a 2 -u\)m 2 . 
Solving for m c , the tangential wavenumber at cutoff, 

v 0 ± Ja z -U 0 


m c = m = © 


2 2 2 
<r-u 0 -v 0 


or 


m c = 


© 


fA/„±/r 


Mi 


1 -M 2 


(A-20) 


Thus the £ and q components of the phase velocity are 
equal to the mean flow velocity components u o and vo. 
This solution corresponds to convection of vorticity by the 
mean flow with no associated pressure fluctuations. 
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Appendix B Frequency relations and 
scattering effects 

Let a disturbance be described as 


to 2 = -mU\ =-§ftA (B-9) 

Waves created by interaction of the gust with the stator must 
satisfy the interblade phase angle imposed by the gust, i.e., 

m q S2 = P 2 +27i<7, <7 = ±1,±2, ... (B-10) 


B = i?exp[/(cD|f + /^ +mrh)] (B-l) 

in a reference frame moving at velocity 

v$, = 0, v^j = U\ with respect to a fixed £,r| coordinate 
system. Let a second reference frame with coordinates 
t] 2 have velocity v^ 2 = 0, v n2 = C /2 with respect to the 
£,t| coordinates. Thus 

(B-2) 

and 

ri - ry + C/,t (B-3) 

for 7=1,2. Using Equation B-3 twice, first for j- 1 and then 
for j- 2, then equating the results gives 


where m q is the tangential wavenumber of a scattered wave 
and 

J 3 2 = ms 2 = 27c| 2.. (B-ll) 

Substituting B-l 1 into B-10 and rearranging yields 

'” , = 2 "(yr + «) (B - 12 > 

Using Equations B-8 and B-l 2, any of these scattered waves 
incident on the rotor will be at frequencies (in the rotor 
frame) 

co* = ©2 +m q U\ 

= (m q -m)U\ 


Til =1l2+(f/2- U\)t. (B-4) 

Substituting B-2 and B-4 into B-l yields 

i? = £exp {/[©]/ + ^2 +/n(ri 2 + (^2 - C/ 1 ) 0 ]} 

= B exp {/[(©j +m(U 2 - Uj))r-f /^ 2 + ^t] 2 ]}. (B-5) 

B-6 implies that © 2 > the frequency in the second reference 
frame, is 

©2 =©1 +m( U 2 - C/j ) . (B-6) 

Thus the frequencies generally differ due to relative 
tangential motion between blade rows, with the frequency 
shift depending upon the tangential wavenumber m in 
addition to the velocity difference. 

The appearance of m in equation B-6 has important 
consequences when scattering effects are considered. As an 
example, consider the first harmonic forced response of a 
stator (row 2) by the wakes from an upstream rotor (row 1). 
For the first harmonic of the wake, the tangential 
wavenumber is 



(B-7) 


= (B-l 3) 

The rotor will scatter the incident waves to wavenumbers 

m* r = m * + %r r = 0,±l,±2, ... (B-l 4) 

^ 1 

= 2 ”te+iQ 


at frequencies according to Equation B-13. Finally, the 
m qj waves coming from the rotor will excite the stator at 
frequencies 

©2 = ©* - m q ' r U\ 


2ti( 1 +r) 
s 1 


U 1. 


(B-15) 


The point of this example is that scattering at the rotor 
affects the frequencies seen by the stator and vice versa. As 
shown by Equation B-15, waves scattered at the rotor (i.e., 
r * 0) return to the stator at a frequency which differs from 
© 2 , the frequency of interest, as given by Equation B-9. 
Upon further scattering, some of wave energy will return to 
the frequency ©2 - this effect is neglected in this paper. 


where s\ is the rotor blade spacing. Since U 2 K), the 
excitation frequency of the stator is, from Equation B-6, 


©2 = ©1 -mU \ . 


(B-8) 


But the wake disturbance frequency is zero in the rotor 
frame so ©1 =0 and 
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Figure 1 Wave propagation from isolated blade row 
due to forced response 


Figure 2 Wave propagation due to forced response of blade 
row in multistage turbomachine 



Figure 3 Cascade nomenclature 
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Figure 4 Rotor-stator interaction noise 
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Figure 5 Configuration I isolated rotor stability (interblade 
phase angle adjacent to symbol) 


Figure 6 Configuration I rotor stability: comparison of 
unsteady aerodynamic coupling (AC) and 
no coupling (NC) solutions - 
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Figure 7 Configuration I stability, dynamic coupling 
(DC) solution 
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Figure 8 Configuration I stability: comparison of unsteady 
aerodynamic coupling (AC) solution with dynamic 
coupling (DC) solution 
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Figure 9 Effect of axial gap on least stable solution, 
configuration I 


Figure 10 Configuration II stator stability: unsteady 
aerodynamic coupling (AC) solution vs. 
no coupling (NC) solution 



Figure 11 Configuration II stability: dynamic coupling solution 
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